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ABSTRACT 

We construct a steady analytic accretion flow model for a finite rotating gas cloud 
that accretes material to a central gravitational object. The pressure gradients of the 
flow are considered to be negligible and so, the flow is ballistic. We also assume a 
steady flow and consider the particles at the boundary of the spherical cloud to be 
rotating as a rigid body, with a fixed amount of inwards radial velocity. This represent s 
a generalisation to the traditional infinite gas cloud model described by lUlrichl (|1976l ). 
We show that the stream lines and density profiles obtained deviate largely from the 
ones calculated by lUlrichl . The extra freedom in the choice of the parameters on the 
model can naturally account for the study of protostars formed in dense clusters by 
triggered mechanisms, where a wide variety of external physical mechanisms determine 
the boundary conditions. Also, as expected, the model predicts the formation of an 
equatorial ac cretion disc ab out the central object with a radius different from the one 
calculated by lUlrichl (| 19761 ). 
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1 INTRODUCTION 

The first model of a spherical symmet ric accretion flow to- 
, wards a central object was made by iBondil (|l952l) . With 
time, it has become the key to under stand different accr e- 
' tion processes in the universe (see e.g. lFrank et~ai] (|2002l )). 
despite the fact that it was created for mathematical curios- 
ity, rather than for astrophysical applications (|Bondill2005l ). 

IBondil thought of an infinite gas cloud with a non- 
relativistic gravitational central object located at the origin 
of coordinates. He assumed the flow to have reached a steady 
state. Under these assumptions iBondl was able to integrate 
the hydrodynamic equations for this spherical symmetric 

flow. 

Years later, lUlrichl (|l976l ) modified iBondH 's ideas by as- 
suming all fluid particles at the border of the cloud to have 
a certain amount of angular momentum, following a dis- 
tribution give n by a rigid body rotating about the z axis. 
In his model, lUlrichl took no account of pressure effects on 
the infalling gas. In other words, his analysis was ballistic, 
which is approximately true if the flow is supersonic and if 
heating by radiation and viscosity effects are negligible (cf. 
lMendoza'et~aIll2004 ). 

The boundary conditions in lUlrichT s model imply that 
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all fluid particles move in parabolas whose focus correspond 
to the origin of coordinates. As soon as particles arrive to 
the equatorial plane they collide with their symmetric coun- 
terparts and stay on it, forming an equatorial disc of radius 
r u ma de of flu id particles which rotate about the central 
object. lUlrichl 's model is usually taken as the basic model 
which predicts the formation of an accretion disc of parti- 
cles orbiting a central object due to accretion with rotation. 
The solution given by lUlrichl \l97^_ was extended by Casen 
fe Moosman and is used (see e.g. lNage]||2007l ; lLin fc Pringlel 
Il99d ; IStahler et al.lfl994l ) as a natural initial condition for 
simulations dealing with the problem of disc formation and 
subsequent evolution. 



Both models, IBondil 1 19521 ) and lUlrichl (|l97rj ) have 
been taken to the relativistic regime by assumi ng the 
centr al object to be a Schwarz schild black hole (|Michell 
1 19721 ; iHuerta fc Mendozal 120071 ). They have wide ap- 
plications in high energy phenomena (see for exampl e 
iLee fc Ramirez-Ruizir2006l ; lBeloborodov fc Illarionovll200ll ). 

On many astrophysical situations where Ulrich's model 
is used, the usual assumption is that the boundary condi- 
tions for an infinite gas cloud remain valid even when it is 
applied to the inner core of a large accretion cloud (see e.g. 



Cassenfc Moosman|[l98ll ; lLin fc Pringl c 1990: St ahler et al.l 



19941 ). This is not necessarily true, because the accretion 
flow at the boundary of the inner core has, in general, an 
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inward radial velocity different from zero3- Also, if the an- 
gular momentum of particles follows rigid body rotation at 
this finite radius then, as we describe in this article, the flow 
may largely deviate from Ulrich's solutions, since the spe- 
cific total energy (kinetic plus gravitational) is not zero and 
so, orbits are no longer parabolic. The relaxation of Ulrich's 
boundary conditions leads to a more realistic infall model 
with applications to a wider variety of accretion problems. 
With an appropriate selection of boundary conditions this 
could in principle modify the main results of a number of 
previous works in the literature dealing with various astro- 
physical accretion problems. 

In this article we assume a central Newtonian potential 
produced by the central star, and a rotating gas cloud with 
a finite radius ro- Based on this, we construct a cylindrical 
symmetric flow which also predicts the existence of an accre- 
tion disc about the central object. As it will be shown, the 
trajectories are no longer sections of parabolas. They are 
general conical sections and once boundary conditions are 
fixed they can be of different families for each fluid particle. 

In order to show the relevance of having a finite cloud 
radius ro and a non-zero radial velocity v ro , note that this 
boundary conditions can be added to the standard Ulrich's 
ones which are defined by the mass M of the central star, 
an accretion rate M and a specific angular momentum ho 
of the cloud. In other wo rds, r p and v ro characterise the 
modifications to Ulrich's (|l97rj ) model. The length ro can 
be taken as the outer boundary of the inner region of the 
cloud, where its self-gravity is neglected, and the main force 
is given by the gravity of the star. 

In real astrophysical problems this model can be ap- 
plied. For example to NGC 1333 IRAS 4 region where 
an unus ual velocity profil e in an inf ailing envelope is pro- 
posed by I Choi et al.l (|2004I ) to explain the observation. Also, 
it is relevant in the study of triggering star formation, 
because whatever the physical mechanism responsible, a 
whole set of combinations for ro and v ro is possible. In 
this scenario, Ulrich's model is not applicable. There are 
well known mech anism for the t riggering, namely stellar 
ejecta from w i nds (Foster fc B oss 1996) or supernova shells 
(Melioli et all l2006n ~ as well as the ionising radiation from 
a star, or a collection of s t ars, that produces an ionisation 
shock front (|Bertoldill 19891 ; lEsquivel fc Ragall2007f l . Each of 
these mechanisms sets particular boundary conditions for 
the collapse of the cloud to form a star. 

The triggering is relevant in clouds with a high density 
of potential in sites of star formation. A first generation 
of stars can influence the dynamical configuration for the 
next generation. This mechanism is highly important when 
it is realized that the stars are formed in turbulent clouds 
l|Vazouez-Semadeni et al.l 120071 ) and that a qui escent core 
for star formation envisioned bv lShu et al.l (jl987h is in some 
dense clusters difficult to fi nd. A complete p i cture of these 
ideas is clearly presented bv lHartmann et alj (|200ll ). 

A glimps e of the full range of bound ary conditions is 
exemplified bv lHartmann fc Burkertl (|2007t ). They speculate 
that the rich Orion nebula cluster is formed by the collapse 
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of an elliptical rotating sheet of gas with a density gradi- 
ent along the major axis. A high density region is formed 
when one tip of the cloud coll apses preferentially and so, 
the already formed protostars (|Heitsch et alj[2008h accrete 
material with velocities given by the global collapse of the 
cloud, independent of the gravity of the collapsed object. 
Thus, free choosing of boundary conditions is highly impor- 
tant to study the regions close to the star. The resolution of 
these regions is beginning to be reached by computer simu- 
lations and observations. 

In this context, the ballistic solution given in this paper 
is useful to study star and disc formation in very complex 
environments, where one has the need to choose boundary 
conditions freely and easily. However, an obvious weakness 
in the model is that at ro, the cloud is rigidly rotating. Thus, 
in the scales where this solution is valid, one either needs to 
invoke a magnetic field a t r > rp that is responsible to en - 
force rigid body rotation (|Mouschovias fc Paleologoulll979h . 
or the fact that the shell feeding the star-disc system ro- 
tates rigidly during some reasonable period of time, in or- 
der to consider this process as quasi-s tationary. However , 
there are many works in the literature dKenvon et al.|[l993l ; 
I Javawardhana et al.l l200ll ; IWhitnev et aL I l2003h that make 
such an assumption. 

In Sections [2] and [3] we calculate the velocity and den- 
sity fields of the flo w res pectively. At the end of Section [3] 
we recover luirichl 's (|l976l ) model. SectionUdescribes the be- 
haviour of the flow for its different parameters and discusses 
its limitations. Finally, Section [5] deals with realistic astro- 
physical situations and the validity of our model on those 



scenarios. 



2 THE MODEL AND ITS VELOCITY FIELD 

Let us assume that a spherical cloud with radius rp accretes 
matter in a steady way towards a central object located at 
the origin of coordinates. All fluid particles located at rp are 
taken to follow rigid body motion in the azimuthal direc- 
tion, i.e. they rotate about the z axis of coordinates. At this 
position, particles also have a radial velocity component v ro . 
Here and in what follows, due to the symmetry of the prob- 
lem we use spherical coordinates r, 6, and <f> for the radial 
coordinate, and the polar and azimuthal angles respectively. 
The particle number density n at ro is taken to be constant 
with a value no- To simplify the problem, we assume the 
self gravity of the gas cloud to be small as compared to 
the gravitational potential produced by the central object. 
In addition, if the gradients of pressure are small compared 
to the kinetic and potential energies, then all fluid particle 
trajectories can be approximated as ballistic. The total an- 
gular momentum of the gas cloud points in the direction of 
the z axis, which contains the central object. Note that when 
ro — > oo and v ro — > , the model c onverges to the accretion 
pr oblem proposed b vlUlrichl l| 19761) . which is well described 
by iMendoza et all ( 2004 1. In fact, Ulrich's convergence is 
mathematically proved at the end of Section [31 

Since the gas is being accreted from ro, it's accretion 
flow rate per unit mass N at any fixed radial distance r is 
constant and given by 



N — 4-Kr nov r 



const. 



(1) 
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At ro , the distribution of specific angular momentum h 
is given by h — ho sin 9o , where 9o is the initial polar angle 
of the particle's trajectory. 

Under all the above assumptions, the specific energy 
E and the specific angular momentum h are constants of 
motion along each particular trajectory. The specific energy 
E is given by 



E 



1 2 I/* 2 

2 Vr + 2^ 



GM 
r 



1 2 I fe 2 sin 2 6> 
2 Vro+ 2 r 2 



GM 

ro 



(2) 



where G is Newton's gravitational constant. We now intro- 
duce two dimensionless parameters a and v given by 



2 

a := 



hi 
r 2 E 



Eo 



(3) 



where r u = ho /GM is the disc's radius in Ulrich's model and 
i?o := GM/r u is the specific gravitational potential energy 
of a fluid particle evaluated at r u . The quantity fi represents 
the ratio of Ulrich's disc radius to the original cloud's ra- 
dius. So, for example, in order to consider an infinite cloud's 
radius we must take the limit u — > 0. On the other hand, 
v represents the initial amount of radial velocity measured 
in units of the Keplerian velocity Vk := VEo, at position 
r u . With these parameters, equation ((2} can be written into 
dimensionless form as 



v 2 + fj 2 sin 2 0o — 2/i, 



(4) 



where the dimensionless energy s is given by e := 2E/Eo- 

The trajectory of each fluid particle is contained on a 
plane and is given by a conic section. The origin represents 
one of the foci of the orbit. At a specific initial position, 
the particle is located at ro and its polar and azimuthal 
angles are given by Go and cj>o- Over the orbit plane, the 
particle trajectory is defined by an azimuthal angle ip which 
at the initial position has the value po- The trajectory of a 
given particle is defined b y the solution to Kepler's problem 
(Landau fc Lifshitz|[l989l ): 



Using equation (j8j) to rewrite equation ((5} yields 



sm Uq 
1 — e cos f ' 



where 



£ = cos 1 



+ <po- 



(9) 



(10) 



v COS #0 

We now use the fact that in spherical coordinates 
v,p — r sin 8d(j>/dt — ft sin #o A sin (9, ve = rd9/dt = 
{d9/d(j}) (tty/sin0), and v r = dr/dt = {dr/dO) (v e /r). With 
this and using equation ([5| the expressions for the velocity 
field, measured in units of axe given by: 



«0 



sin 2 6o 
r sin 6 ' 



sin Qq , 2 
vg = — : — - (cos t>o — cos 6M 



2 MV2 



e sin ^ sin 9o 
r(l — ecos^) 



(11) 



(12) 



(13) 



3 PARTICLE NUMBER DENSITY 

In order to obtain the particle number density field we start 
from the continuity equation for a stationary flow, i.e. 



V • (nv) = 0. 



(14) 



We now integrate the previous equation over a volume 
consisting on a collection of streamlines. In other words, this 
"tube" is built in such a way that its lateral surface is bound 
by streamlines and its lids by spherical sections, the upper 
one at the cloud's border, i.e. r = ro = 1/fi and the lower at 
any arbitrary radial distance r such that < r < 1/fJ,. With 
this selection, equation (|I4p can be integrated to obtain: 



sin fo 

— e cos ip 



(5) 



In the previous equation and in what follows, unless stated 
otherwise, distances are measured in units of r u . The eccen- 
tricity e of the orbit is given by 



e = \/ 1 + esin 2 (fa- 



(6) 



At the border of the cloud r = ro = 1//U and so, substitution 
of this in equation ([5]) gives the following condition for tpo m . 



COS Ifio 



— (1 — fis'm 2 do)- 
e 



(7) 



Performing standard spatial rotations its easy 
to obtain the following formulae between the angles 
ip, -po, 9, Go, 4> and <j> : 



nv ■ da 



— nv ■ da 



(15) 



r=l/fi 



Taking into account the fact that the differential area ele- 
ment da is given by 

da = r 2 sin 6 d9 d<f> e r + r sin 9 dr d(f> eg + r dr d9 e^, (16) 
then, eauation (|15[) transforms to 



sin So d#o d(f)o = —nv r r sin# d9 d(j>, 



(17) 



where the particle number density n is measured in units of 
n u := N /4-KV^r 2 . Using the symmetry of the quantities in 
the azimuthal angle involved in the previous equation, we 
can integrate equation (|f 7[) with respect to <j> from to 2n, 
to obtain 



cos(v3 — ipo) = 



cos 9 
cos 9o ' 



M = 



tan 6 
tan ( 



(8) 



I sin So 
r 2 sinS 



-V r 



88 
d9o~ 



(18) 
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On the other hand, from equation (|9} it follows that 



sin(9 / 80 



sin Oq 



sin Oq 



1 + 



2/x cos Oq cos 9 ) I 1 



3 cos 2 Oq — 
r 

2 a \ V 2 



COS 2 #0 



+ 



(19) 



+ 



Sin Oq 



fl + cos 2 9- 2 cos 2 9 ) 



Substitution of equations (|19|) and (|13[) into (|18[) leads to 
the required particle number density field: 



sin^o f /, 3 cos 9a — 1 
n = — — <! I 1 H 



2/i cos 9q cos 1 I 1 — 



cos 2 8q 



+ -r^r ( 1 + cos 2 6»-2 cos 2 9 ) 



1/2 



+ 



(20) 



sin 6*o 



As a result of the collision at the equatorial plane be- 
tween the streamlines coming from the northern hemisphere 
with those coming from the southern one, a shock front is 
expected to form, as well as t he formation of an equatorial 
disc shape structure (see e.g. lUlrichl Il976l ; iMendoza et al.l 
|2004| . and references therein). 

Taking 9 — > n/2 in equation (|9]) followed by the substi- 
tution Oq — > 7r/2 we obtain an expression for the accretion 
disc radius given by: 



1 



(21) 



From this last equation and since naturally < ro, the 
following condition between the parameters arises: 



H < 1 + v 



(22) 



If in equations @, jnj, |T2")). (O and JgD} we take the 
parameters values /x = and v = , we get the following 
set of equations 



sin Oq 



1 — (cos 9 1 cos Oq) ' 



sin 9q ( cos 00 + cos f 



COS Oq — COS f 



1/2 



sin 9 , 2 a 2 mi/2 

t'e = — : — t (cos 6>o — cos 9) . 



V4, 



sin 2 6*o 
rsin6> ' 



(23) 



(24) 



(25) 



(26) 



sin 9q f_ 3 cos Oq — l\ /, cos" 
n= — 1H 1 



COS 2 0q 



-1/2 



, (27) 



which corresponds to the expressions calculated by lUlrichl 
ll 19761) . 




Figure 1. The figure shows projected streamlines of the flow 
for different values of the parameter v (see text) at an arbitrary 
azimuthal angle <f> = const in such a way that the parameter 
H = 1.28. From left to right, top to bottom, these values are 
v = 0.28, 0.8, 3.2, 12.8. These projections do not correspond 
to the real streamlines of the flow (see section |4j . Distances are 
measured in units of the cloud's border r = tq = 1/fi. 



4 STREAMLINES AND DENSITY PROFILES 

Figure [T] shows projected streamlines for different values of 
the parameters fi and v with height z and radial coordinate 
R := rsin^. Each dot on a fixed curve correspond to the 
intersection of the real streamline and the plane (j> — const 
as it is being swept from cj> = to c/> = 7r/2. In other words, 
each line is the intersection of a family of real streamlines all 
having the same Oo, with a fixed plane cj> — const. The rea- 
son of doing this is because direct projections of streamlines 
on a plane <j> = const show false intersections of the stream- 
lines on the obtained plots. Figure Q] shows the projected 
streamlines for a fixed fi and different values of v. 

Figure [2] shows also projected streamlines of the flow 
for a disc's radius of one half of the original cloud's radius 
(t& = 1/2/i). The bottom panels in Figure [2] have "trian- 
gles" drawn on them. These are zones out-of-range in the 
model due to real intersections of the accretion streamlines. 
Note that this behaviour appears because the model is bal- 
listic. We have explored numerically in which cases these 
intersections occur and it turns out that for r& = 1/2^, the 
streamlines show intersections for v > 1. From values of 
v > 50 intersections appear but the out-of-range region is 
kept fixed. 

Figures [3] and [4] show density isocontours and density 
profiles for a fixed height measured from the equatorial plane 
respectively. The out-of-range region appears only in the 
bottom panels of Figure|3] although on the left-bottom panel 
it appears as a tiny region. 

Figure [S] shows plots of the streamlines for the case 
v — 0. The streamlines on the plot have different values 
of the initial cloud's radius. It is interesting to note that the 
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Figure 2. The figure shows projected streamlines of the flow for 
different values of the parameters fi and v (see text) at an arbi- 
trary azimuthal angle = const in such a way that the radius of 
the disc has a value rj = 1/2^. From left to right, top to bottom, 
these values are p, = 0.5, 0.81, 1.31, 51.01, v = 0, 0.64, 1.62, 101. 
These projections do not correspond to the real streamlines of 
the flow (see sectional. The filled "triangles" drawn in the bot- 
tom plots are zones out-of-range where streamlines of the flow 
intersect. This occurs because the flow has been assumed to be 
ballistic. Distances are measured in units of the cloud's radius 
r = ro = 1/fJ.. 
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Figure 3. The figure shows logarithmic particle number 
density n isocontours in units of n u . From left to right, top 
to bottom the values of the parameters of the flow are the 
same as in Figure [2] and the values of the isocontours for 
each panel are given by {0.32, 0.4, 0.6, 0.75, 1.5, 1.75}, 
{0.3, 0.34, 0.4, 0.5, 0.65, 0.9, 1.3, 2.0}, 
{0.0, 0.13, 0.33, 0.45, 0.6, 1.0, 1.5, 2.0} and 
{-1.0, -0.85, -0.65, -0.45, -0.32, -0.2, 0.0, 0.5} re- 
spectively. The dark regions on the bottom panels are the 
out-of-range zones. Distances are measured in units of the 
cloud's radius r = rg = 1/fi. 



different streamlines having the same Go as initial condition, 
arrive to the same projected position over the disc. 

Figures [6]|8] are equivalent plots as those of Figures 12141 
for the case /j, — (i.e. ro — > oo), which c an be thought 
of as the first natural modification to lUlrichl 's model, all of 
them depending on the single dimensionless parameter v. 
Note that in all these cases, there is no intersection of the 
flow streamlines. This appears to indicate that as long as 
the cloud reaches a sufficiently large radius ro, with respect 
to the accretion disc radius r<i, the intersections disappear. 
Indeed, we have tested different accretion models varying 
Td/fo from ~ 0.0 to 1.0 and the previous statement turns 
out to be correct. As an example, a value of r^/ro ~ 0.08 
represents the threshold of intersections of streamlines for 
values of v = 10.0. 

From Figures[U]and[7]it is clear that the particle number 
density does not accumulate that much in the origin as v 
increases. 



5 ASTROPHYSICAL APPLICATIONS 

From the point of view of Spectral Energy Distribution 
(SED) modelling of protostellar cores, one of the main in- 
gredients required for a radiative transfer calculation is 
a specification of the density profile in the core. Many 
authors use Ulrich 's model to achieve this (see for e.g . 
iKenvon et alll 19931 ; iButner et al1ll99ll ; I Adams fc Shulll98rj ; 



IWhitnevet alll2003r i. A radiation transfer simulation re- 
quires a lot of information: luminosity of the star, compo- 
sition, abundance and optical properties of the material in 
the core, in terms of its location with respect to the star. 
Thus, a full analysis of the SED produced by the use of the 
density given in this article is beyond its scope, but must be 
consider in the future. 

In order to give a qualitative idea of the difference in 
the results, using either Ulrich's density or the density cal- 
culated in Section [3] of this article, we note that less density 
means more luminosity arriving from the object. The main 
argument is that the optical depth, which is proportional 
to the surface density calculated including all the material 
along a line of sight, is an element that give us an idea of 
the SED produced by the core. 

For a disc viewed pole-on, in Figure [9] we present the 
axisymmetric surface density £ as a function of the radial 
coordinate, for Ulrich's model (E u ) and from this paper (E a 
for v ro = and Sj, for v rn = 5 x 10 4 cm/s, according to the 
results o f lHennebelle et al.ll|2004h ) . This comparison is done 
for 3 representative cases of ro and fl = ho/ro, such that 
r u = 20 AU is the same for all. We take three different sets 
of boundary conditions: (a) ro = 300 AU, Q = 10~ 11 s _1 , 
(b) r = 3000 AU, n = HT 13 s -1 , and (c) r = 30000 AU, 
S7 = 10~ 15 s _1 . Note that in all three cases, ho is kept fixed 
since for all models the radius r u has been chosen to have a 
fixed value. The other common parameters on these models 
are given by the mass of the central star M = 1 M© and the 
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Figure 4. Both panels show density profiles for a fixed azimuthal 
angle, i.e. <j> = const and a fixed height z = 0, 0.1, 0.2, 0.3, 0.5 
measured from the equatorial plane. The figures were drawn using 
values of v = 0.64, 1.62 respectively for a fixed value of the radius 
of the disc = l/2/i. Note that for a height 2 = the particle 
number density diverges at r = 0, ±0.5. Distances are measured 
in units of the cloud's radius r = 1/fA. 
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R 

Figure 5. The plot shows a comparison between the streamlines 
for three different values of the initial cloud's radius in the case 
v = 0. These values are ro/r u = oo, 2.46, 1.38 for the cases of 
continuous, dotted and dotted-blank lines respectively. Distances 
are measured in units of r<j . Bear in mind that for different values 
of tq the disc radius r d varies and so, if one is to keep r d = 1 
then different values of the mass of the central object and the 
angular momentum of the cloud need to be properly readjusted. 
Distances in the plot are measured in units of the radius of the 
disc 7\j. 




Figure 6. The plots show projections of streamlines (in the same 
sense as in Figure [T} which correspond to an infinite value of 
the cloud's radius rrj. From top to bottom the parameter v has 
values of 0, 2, 100 re spectively (see Section [2}. The top panel 
corresponds to lUlrichl 's accretion model. Distances are measured 
in units of the radius of the disc r d . 

mass accretion rate M — 10~ 6 Mq yr^ 1 . The most typical 
as tronomical case is (b) in accordance w i th th e observations 
of iBenson fc Mversl (|l98Sh : Ijiiina et al.l (|l999h . 

In Figure[9l one can see that Ulrich's surface density E u 
is smaller than the one calculated using our model, in the 
inner regions of the core (few times rd). The difference be- 
tween both curves increases as r$ decreases. Thus, the flux of 
radiation in the typical frequencies produced in this region 
are larger in the former than in the latter. Of course, a full 
analysis of this problem awaits an appropriate simulation, 
but since case (b) is quite common in accordance to obser- 
vations, we can see that the modification to Ulrich's model 
is quite important. From the figures it follows that even if 
v ro has a null value, our model shows significant deviations 
as compared to the ones predicted by an Ulrich approach. 

The last part of this section is intended to show 
the relevance of having flexible boundary conditions 
for the accretion flow studied in this article. The col- 
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Figure 7. The panels show the logarithmic particle number 
density isocontours measured in units of n u corresponding 
to the same set of parameters chosen in Figure [6] Dis- 
tances are measured in units of the radius of the disc r<j 
and from top to bottom the sets of contours take the fol- 
lowing values: {-0.55, -0.47, -0.35, -0.2, 0.0, 0.5, 1.0}, 
{0.25, 0.4, 0.6, 0.8, 1.1, 1.5, 2}, and 
{1.7, 1.8, 1.95, 2.1, 2.25, 2.5, 3, 3.5}. 



lapse of non-ro t ating clouds w as st udied by iLarsonl 



Jl969h: iPenstonl dl969f): iHunterl (Il977i ). As an example, 
IWhitworth fc Ward- Thompson! (|200ll ) built a model, in or- 
der to fit lifetimes and accretion rates for class protostars 
(this stage ends until half the final mass of the central object 
is assembled). In this model, the density diff ers from the sin- 
gular isothermal sphere (SIS), cf. IShul ( 19771 ). because in the 
former, the density has an inn er flat profile. Moreover, obser - 
vations of protostellar cores (|Ward-Thompson et ail 119991 ) 
suggest that initial conditions for protostellar collapse de- 
part from SIS. 

At this point, it is clear that the boundary conditions 
in a collapsing core depend on the environment. Dynam- 
ical collapses are present in regions like p Ophiuchi and 
Perseus clusters, where the collapse and subsequent star 
form at ion is triggered b y ext ernal agents (|Motte fc Andrei 
l200ll ). iHenriksen et al.l (|l997l ) stress that strong accretion 
observed in p Oph is related to the impact of an interstellar 




Figure 8. From top to bottom, the figures show the logarithm 
of the particle number density of the flow in units of n u , cor- 
responding to the same values of v that Figure \E\ has. All plots 
have density as a function of the equatorial radial coordinate R 
for different fixed heights z = 0, 0.15, 0.3, 0.5. In all figures, the 
plot at z = has density divergences at R = 0, ±1. 



shock wave. In regions like Taurus it is possible that cores 
collapse by themselves without an external perturbation and 
so, the mas s accretion rate is closer to the value calculated 
for the SIS (|Shulll977l ). 

Triggered star for mation by an exte r nal m echanism is 
throughfully studied bv lHennebelle et all l|2004h . This work 
presents numerical simulations of a rotating core, with high 
spatial resolution, for the study of the formation of a disc. 
The relevance of the perturbation is measured with their pa- 
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Figure 9. The figure shows axisymmetric surface density £ plots 
as a function of the radial coordinate R for a disc viewed pole- 
on, using an Ulrich's model for which the surface density is 
given by S u and from the results of this article with S a for 
cases with v ro = and for v rg = 5 X 10 4 cm/s. Prom top 
to bottom, we have chosen three examples in such a way that 
r u = 20 AU and the corresponding sets of boundary conditions 
are given by: (a) r Q = 300 AU, Q = lO" 11 s -1 , (b) r = 3000 AU, 
Q = 10 _13 s _1 , and (c) r = 30000 AU, Q = 10~ 15 s' 1 . The 
other common parameters on all models are given by M = 1 Mq 
and M = 10 -6 Mq yr _1 . Distances on the plot are measured in 
astronomical units and the surface density in g/cm~ 2 . 



rameter cj>, which is the ratio between the time-scale in which 
the internal pressure doubles, to the sound-crossing time of 
the core. A high value of cj>, represents a slow compression 
(spontaneous collapse); the opposite is found with a small 
value, i.e. fast compression (strongly externally triggered). 

The difference between both cases is clearly seen in the 
size-scale of the disc. In the first one, the disc evolves to a 
two-armed spiral pattern without fragmentation. Ring for- 
mation is the outcome usually found in the second case. 
Most of the material reaching the centre has too much 
angular momentum and so, it m oves outwards and estab- 
lishes a dense ring l|Nagell l2007h . This configuration pro- 
motes fragmenta t ion. D espite the fact that in this case 
IHennebelle et al.l (|2004 ) found no central protostar, there 
must be situations where certainly it is present. In conclu- 
sion, either spontaneous or externally triggered collapse can 
be studied with the ballistic s olutio n presented in this paper. 

Both, IHennebelle et al.l (12004 ) and our work start from 
a core in solid body rotation. The former evolves to a core 
of differential rotation, where a rapid compression implies 
higher velocities and densities in the outer parts of the disc. 
This fast compression intuitively means that a larger v rQ 
is present. Indeed, if in our model we set a large v ro , then 
higher velocities and densities are obtained in the same re- 
gion. For, substituting 9 = n/2 and 9o in terms of r from 
equation ((9]) in equation (|I3I) it is found that for a given 
r, the velocity v r and particle number density n increase 
with an increasing v. This conclusion, expected intuitively, 
reinforces the idea that the solution presented in this paper 
can be taken for studies of disc formation with an analytical 
approach. 

The robustness of this kind of solutions is that, depend- 
ing on the particular astrophysical situation, we can choose a 
ballistic solution with appropria te boundary conditions. For 
example. lThroop fc Ballvl (|200&f ) studied numerically the ac- 
cretion of material into a star-disc system moving through 
the gas in a cluster. An appropriate solution, should help to 
study small scales, where the star's gravity represents the 
main cont ribution. For the case o f star formation in molecu- 
lar clouds. iBate fc Bonnelll (|2005h describe the accumulation 
of material as "competitive accretion". The matter comes 
from different places with different boundary conditions. In 
this case, ballistic solutions are able to give glimpses of the 
physics behind the formation of star-disc systems. 



6 CONCLUSION 

We hav e constr ucted an analytic accretion flow that expands 
lUlrichl 's (| 19761 ) model in the sense that the radius of the 
rigid body rotating cloud is finite and an allowance for a 
radial accretion velocity at the cloud's border is assumed. 
When the radius of the cloud tends to infinity and the radial 
velocity of the input flow at infinity goes to zero, the model 
converges to that of lljlrichl . 

We have also shown that for real astronomical systems, 
where typical sizes of the clouds are ~ 30000 AU and for 
which the initial inward velocity is null, the differences with 
Ulrich's model is quite noticeable. The problem with this 
hypothesis is that to achieve a fixed accretion rate, then 
the density at the cloud's border needs to be infinite. The 
difference between both models becomes even greater if one 
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assumes an initial inward velocity different from zero. Since 
this assumption makes the density not to diverge at the 
clouds border, it is the more reasonable model to take for a 
real astronomical system. 

We are developing a relativistic model with all these 
properties in order to gen eralise the model constructed by 
iHuerta fc Mendoza lj2007f ). These results will be the subject 
of a subsequent article. 
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